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Abstract —In order to foster electric vehicle (EV) adoption, 
there is a strong need for designing and developing charging 
stations that can accommodate different customer classes, distin¬ 
guished by their charging preferences, needs, and technologies. 
By growing such charging station networks, the power grid 
becomes more congested and, therefore, controlling of charging 
demands should be carefully aligned with the available resources. 
This paper focuses on an EV charging network equipped with 
different charging technologies and proposes two frameworks. 
In the first framework, appropriate for large networks, the EV 
population is expected to constitute a sizable portion of the light 
duty fleets. This which necessitates controlling the EV charging 
operations to prevent potential grid failures and distribute the 
resources efficiently. This framework leverages pricing dynamics 
in order to control the EV customer request rates and to provide 
a charging service with the best level of quality of service. 
The second framework, on the other hand, is more appropriate 
for smaller networks, in which the objective is to compute 
the minimum amount of resources required to provide certain 
levels of quality of service to each class. The results show that 
the proposed frameworks ensure grid reliability and lead to 
significant savings in capacity pfanning. 


I. Introduction 

A. Motivation 

Electric Vehicles are becoming a viable transportation op¬ 
tion as they offer solutions to an array of current societal 
problems ranging from high oil prices to environmental con¬ 
cerns. Corollary, EV population is expected to reach a sizable 
market portion in the next decade (e.g., 10% of the U.S. 
National fleet and similar targets in Europe) 0. However, 
achieving such penetration rates requires wide deployment of 
charging facilities that can serve different types of charging 
requests. On the other hand, if not controlled, EV charging 
can easily lead to transformer and line overloading. Moreover 
it can deteriorate the power quality and even endangers the 
security of supply. The impacts of EV charging on the grid 
has been very well documented in (TJ— O and 0. For instance 
0 argues that if 5% of the EVs charge simultaneously using 
fast charging technology, 5.5 GW of extra power would be 
needed in Virginia and Carolinas region by 2018. Similarly, 
NERC regions would require an expansion of 5.5% in their 
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power generation capacity in a typical EV penetration rate 
of 25%. Furthermore distribution grid could easily become a 
bottleneck. For example, authors of 0 state that even adding 
two Level-2 chargers in a typical neighborhood in the US 
could easily cause transformer overloading. Also uncontrolled 
demand can decrease the efficiency of the power grid opera¬ 
tions and increase the power generation cost if it occurs during 
peak hours. 

B. Related Works 

There has been an increasing body of literature on de¬ 
veloping intelligent charging station architectures a-im 
controlling and scheduling EV demand m. Ha-ini and 
rather sparse literature on capacity planning on charging 
infrastructures serving multiple classes of customers m, ED. 
This section provides a brief overview of related literature. 

The studies on charging station design can be divided 
into two categories. The first category includes works from 
power engineering perspective 0-0 where the studies aim 
to minimize the charging duration by improving the efficiency 
of power electronics and aiding the system with energy storage 
system. On the other hand, approaches in the second category 
are mostly focused on the system level where they abstract 
the details of the underlying power system components and 
evaluate the system performance in terms of long term sta¬ 
tistical metrics e.g., mean waiting time and percentage of 
served customer using the arguments from queuing theory 0, 
0, 02. Study in 0 presents a small scale fast charging 
station design and blocking probability, that is the probability 
that the request of an EV is rejected, is used as the main 
performance metric to solve the optimal resource provisioning 
problem. In iflOl authors use M/M/s queues to model the EV 
demand at fast charging stations near highway exits. In 0 
residential charging infrastructures are modeled using M/M/oo 
queue, where they broadcast the arrival rates to customers 
such that probability of blocking (exceeding circuit capacity) 
is zero. Furthermore, the work in 02 uses BCMP queueing 
network model to estimate the EV charging demand. In this 
paper we perform a more holistic approach by taking into 
account that the charging station can serve multiple classes 
of customer that are differentiated by the charger technology 
used. Also, we employ loss-of-load-probability (LoLP) as our 
main performance metric, which measures the probability that 
grid resources cannot accommodate EV demand 11201 . 

The pricing-based control methods have been employed 
successfully in communication networks where the goal is 
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Fig. 1: Two design problems. 


to match scarce bandwidth resources to insatiable user de¬ 
mand ED. In a similar fashion, pricing mechanisms appli¬ 
cations gained popularity in the literature on controlling EV 
demand m, m-m and l22l . Authors in (T) and ED 
propose a control mechanism in a network of charging stations 
to route customers to idle stations. The work in lfl4l proposes a 
pricing scheme for EV chargings that leads to socially optimal 
solution, whereas ll22ll uses proportional fairness pricing from 
communication literature to control the charging rates of EVs. 
In our work we also consider the technological constraints 
and current state of affairs that charging rate is constant 
until the end of the service. The work presented in M 
is most relevant to our work, in which researchers model 
general power consumers behavior with utility functions and 
propose a pricing algorithm to control the consumption of 
smart grid users. Our work mainly focuses on multi-class EV 
chargings and our framework ensures grid reliability at all 
times. The literature on resource provisioning is rather sparse. 
The study in lfl8l presents a capacity planning framework in 
a large scale charging stations with single class customers. 
This framework is based on computing a deterministic quantity 
“effective power” using On-Off Markov models. In this study 
we assume that charging infrastructure can serve multiple 
types of customers. 

The proposed methodology is rooted in multi dimensional 
loss systems in teletraffic engineering, where the goal is to pro¬ 
vide statistical quality of service(QoS) guarantees to customers 
with different demand profiles. The control mechanisms with 
congestion pricing in multi-rate Erlang-B systems and related 
resource provisioning problems are addressed in l23l - |Z7| . 
The works presented in f26\ and l22l . J24|, lf25ll focus on 
efficient computation of blocking probabilities and the their 
derivatives, whereas studies in |23l and j27| provide efficient 
algorithms to solve bandwidth provisioning problem in con- 
gestible networks. 

II. Proposed Frameworks 

Electric vehicle charging demands are primarily dominated 
by customers’ needs and preferences. Based on this premise, 
designing different aspects of the network of charging stations 
is governed by the economic dynamics between network 
operators, on one hand, and the customers on the other hand. 
In such dynamics the former seek operational reliability and 


profit margins and the latter seek economic incentives. Such 
dynamics and their pertinent design challenges vary by the 
size of network and customers population. In this paper we 
consider two canonical design problems for EV charging 
infrastructures fed by a single substation suited for large 
metropolitan areas and small cities as summarized in Fig. Q] 

Specifically, the first framework proposes a pricing-based 
EV control for charging stations located in large metropolitan 
cities in which the EV population is assumed to be large. 
In this framework, the objective is to incentivize and control 
the customers to submit their charging demands at certain 
optimal rates which maximize a social welfare utility. This 
pricing-based approach framework, at its core, enables station 
operators to alleviate the congestion and the degradations in 
the service quality during increased demand periods through 
incentivizing the customers to defer their charging needs 
to less congested periods. This control mechanism ensures 
charging services with desirable level of QoS guarantees. 

The second framework focuses on small cities, in which 
customer demand can be obtained from profiling studies, 
and proposes a capacity provisioning mechanism for EV 
charging stations. In contrary to large metropolitan areas, small 
cities with fewer charging stations and customers enjoy less 
flexibility for incentivizing the customers. Hence, the primary 
goal in this framework is to compute the minimum amount of 
required resources to provide charging service to meet a target 
level of QoS. 

A. Network Model 

We consider a network of charging stations serving N cus¬ 
tomers. Since the customers do not necessarily seek identical 
charging services, we consider J distinct customer classes, 
which are distinguished by their preferences; size of the bat¬ 
tery packs, amounts of requested demands, and the available 
charger technologies. The charging rate of customers of type 
j £ {1,.., J} is denoted by bj. Since the amount of resources 
are finite, upon the arrival of a specific customer type j, 
if the available resources is less than bj the customer will 
not be served resulting in an outage. Hence, the probability 
of being blocked (or loss-of-load-probability) constitutes a 
natural system performance metric. Due to the temporal fluc¬ 
tuations network conditions (e.g., congestion etc.), the amount 
of available resources varies over time. Hence, we consider a 
dynamic system indexed by the time index k £ IN in which 
Ck denotes the aggregate units of grid power available to 
the entire EV fleet at time k. Note that Ck is an exogenous 
parameter given to station. Since the charging stations reside 
in a small well-confined regions, the set of chargers can be 
abstracted collectively as a super-station with multiple classes 
of customers. We further assume that system capacity is large 
when compared to demands of each EV type. 

Increased outage events leads to reduced utility for the 
network operator and service disruption for the customers. 
Hence, optimizing the operation of the network strongly hinges 
on finding the decision rules that lead to the optimal level of 
Loss-of-Load-Probabilities (LoLP) that is defined as the prob¬ 
ability that grid resources fall short of aggregated customer 
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Fig. 3: Customer preferences can be determined by EV type. 
Fig. 2: Pricing-based control mechanism for social welfare amount of reque sted demand, and the charger technology, 
maximization 


demand {20fl . The leverage that the network operators have for 
adjusting the LoLP rates at the desired optimal rates is service 
pricing, by using which they can influence the behavior of the 
customers abstracted by their arrival rates. To formalize this, 
we define Pj(k) as the price of service of type j and define 
the price vector p(k) = [pi(fc),... ,pj(k)], accordingly. 

Based on the prices p(k) and its needs, each customer 
decides whether to generate a service request of a certain type. 
We denote the rate of type j G {1,.., J} requests generated 
by customer n G {l,..,iV} during period k by A"(fc;p(fc)). 
Hence, the aggregate arrival rate of service requests of type j 
is 

N 

A j(k;p{k )) = ^2 Xj(k-,p(k)) . (1) 

n= 1 

From the network operators point of view, these requests 
are being constantly generated across the networks and as 
discussed in m, m, ma, a common assumption is that such 
aggregated requests (and not necessarily individual requests) 
arrive according to a Poisson process. Accordingly, we define 
the arrival rate vectors 

A n (k-,p(k)) 4 [X?(k;p(k)),...,\ n j(k-,p(k))} , (2) 
and X(k;p(k)) = [Ai(fc;p(fe)),..., Xj(k\p(k))] . (3) 

By noting that the FoFPs for each class are functions of 
the arrival rates, finally we define /3j(k; X(k;p(k)) as the 
probability that customers of type j are blocked and define 
P(k;X(k;p(k)) as the corresponding FoFP vector. 

Since different customer types do not necessarily charge 
their batteries at the same rate, we denote the average charging 
duration of customers of type j by 1/pj. It is noteworthy 
that from station operator standpoint, the information required 
to characterize customers’ profiles are the average service 
duration and the charging rates. These parameters enable 
modeling customers’ behaviors with fine granularity. 

The main goal is to keep the aggregated demand below a 
constant power Ck with minimal or controlled loss-of-load 
(outage) events. Drawing constant power is highly desirable 
from the viewpoint of the network operators and its benefits 
include (1) grid components are isolated from stochastic 
variations and hence grid reliability is ensured fT], (8); (2) 
station operator can sign long-term contracts and benefit from 
the lower prices; (3) constant demand will reduce the peak- 
to-average demand ratio of the whole power system and 


accordingly the average spot prices would reduce; and (4) it 
leads to more efficient market equilibrium HI. It is noteworthy 
that since the charging stations reside in a small well-confined 
region, the whole set of chargers act as one big station with 
multiple classes of customers. Due to this fact, power system 
losses are assumed to be negligible and all customers share 
the same resource pool. 

For the simplicity of the notations, in the rest of the paper 
the explicit dependence on k is omitted. 

B. Large Networks: Pricing-based Control Framework 

As the EV population is expected to constitute a sizable por¬ 
tion of the national light duty fleets, there is a need to control 
EV charging operations in order to protect the grid components 
from potential failures. To this end, this framework focuses on 
charging infrastructures that are located in big metropolitan 
areas with high EV penetration rates. The primary goal is to 
leverage pricing in order to control the EV customer arrival 
rates such that station operator can provide a charging service 
with a small level of customer LoLP QoS for each class. Here, 
the primary goal is control the large EV demand in a big city 
for a given amount of grid power (Ck), which is assumed to 
computed by the utility operator and meets the power network 
constraints (e.g., congestion, transformer ratings etc.). Pricing- 
based control scheme is depicted in Fig. [2] 

Specifically, the aim is to design a distributed control 
framework for maximizing an aggregate utility (social welfare 
measure), which requires that the EVs operate at an optimal set 
of arrival rates. In order to establish the tools for the distributed 
design let 

U n ( A";/3(A)) 

for a given set of arrival rates A" and block probabilities (3. 
Hence, the aggregate utility in the network is 

N 

R = Y J U n ( A”;/3(A)). (4) 

n— 1 

Therefore, the social welfare problem can be defined as 
the problem that maximizes the aggregate utility R over all 
possible choices of the arrival rates {A”}^ =1 , i.e., 

N 

max Vt/ n (A";/3(A)) . (5) 

{a\...,a w } ^ 

n =1 

We assume the utility function U n , is increasing in the arrival 
rates {A.,}, and decreasing in the LoLPs j3j. Furthermore, 
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Fig. 4: Numerical Evaluation of multi-class EV LoLP. Each Fig. 5: Numerical evaluation of multi-class EV LoLP for 
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it is assumed that U n is concave in A j and continuously 
differentiable in all of its arguments. This problem is treated 
in Section Hill 


C. Small Networks: Resource Provisioning Framework 


In the resource provisioning problem, we are interested in 
computing the minimum amount of grid resources C such 
that system operator guarantees certain levels of reliability for 
serving different types of users by enforcing B :/ < 5j for each 
customer type. This structure suits the resource provisioning 
problems in small cities in which customer arrival rates can be 
obtained from profiling studies with reliable accuracy. Based 
on the constrains defined for LoLPs, the problem can be cast 
as 


C* 


min C 

s.t t 3j < Sj , for G {1,..., J} 


( 6 ) 


This problem is treated in Section [V] 


D. Toy Example 

We provide a toy example before the case studies in order 
to put some of the details into perspective. Even though the 
objectives of the proposed frameworks are different, computa¬ 
tion of LoLP is carried out for a fixed parameter setting (e.g., 
A, p etc.), the methods for computing blocking probability is 
valid for both cases and details are given in Section [V] Let us 
assume that the charging infrastructure draws C = 1000 units 
of power from the grid and serves three types of EV customers, 
i.e., j G {1,2,3}. Customer classes are differentiated by the 
charging technology they use, and by mimicking the current 
charging standards (fast charging, level-II three and single 
phase). It is assumed that b = {50, 7, 5}, the charging rates 
are p\ = 3, p% = 0.42, and p 3 = 0.2, and the arrival rates for 
each type are Ai = A2 = A3 = 14 (on average 14 customers 
arrive at each time unit). Consequently, the resulting LoLPs 
are /3{ OOQ =0.0152, 0\ ooo = 0.0015, and /3f 0 oo = 0.0011. 
Note that customers of type 1 have the highest LoLP mainly 
because grid resources are shared equally among all customer 
classes and type 1 customers use fast charging technology 
which draws more power compared to other classes. Next, 
we evaluate the same system for a wide range of arrival 


rates. We assume that the sum A = ]Cj=i A' var i es in the 
range [1,80] and A j = A/3 for j G {1,2,3}. The results 
are depicted in Fig. [4] It is important to notice that, for 
the given amount of grid resources (in our case grid allows 
C=1000), stations can provide a very good service for light 
traffic scenarios in which the aggregated arrival rate varies in 
the range A G [1,50], which can be a very typical setting 
for small cities. It is noteworthy that instead of providing 
an extremely stringent LoLP constraints (near zero outage 
probability), network operator can back off from (7=1000 to a 
smaller level C < C while guaranteeing a high level of LoLP 
(e.g. 1% LoLP) for each customer classes. 

On the other hand, in the heavy traffic regimes (e.g., A > 50) 
which can occur in large metropolitan cities, in order to 
provide good QoS, either charging resources ( C ) should be 
increased or the arrival rates should be controlled. To shed 
more light on this, we evaluate the system performance for a 
fixed set of arrival rates (Ai=A 2 =A 3 = 10 ) and for a range of 
serving capacity values C G [500,1000] . As clearly depicted 
in Fig .0 more customers will be accommodated as the amount 
of resources increases. Hence, besides upgrading the grid to 
have increased capacity, which might not be economically 
viable in the short term, an alternative approach is to provide 
charging services based on the available resources and with 
the best possible QoS via controlling customer arrival rates. 
The pricing-based control framework discussed in section [ 6 (c)] 
attacks this problem. 

III. Pricing-Based Control 
A. Global Problem 

This treatment provided in this section has a two-fold 
purpose. Primarily, it provides an optimal LoLP policy in order 
to maintain an optimal social welfare utility as characterized in 
©. Secondly, it proves that the solution provided is amenable 
to distributed implementation. This is of paramount significant 
as in an EV charging network the EVs are autonomous entities 
and operate based on the common information known to all 
the EVs (e.g., charging prices) and their local perception about 
network dynamics. Such distributed nature of the EV networks 
necessitates solving 0 in a distributed manner. 

As it will be shown, an important feature of the proposed 
framework is that each EV adjusts its arrival rates solely based 
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on the updated prices that are dynamically announced by the 
station operators, and such distributed adjustment of arrival 
rates by the EVs leads to an optimal global welfare for the 
entire network. 

By recalling that utility function U n is increasing in the 
arrival rates {A-/}, decreasing in the LoLP f3j, concave in 
A j, and continuously differentiable in all of its arguments, the 
aggregate social welfare R = Yln=i U n is is maximized when 
Vn e {1,..., N} and Vj € {1,..., J} 


d R 
dXf 


dU n 


N J 


v V — 

1=1 s=l rs 



(7) 


Since the LoLP of each charge type depends only on the sum 
of the arrival rates of that type A j = J2n=i • from £7} we 
immediately have Vn £ {1,... ,N} and Vj £ {1,..., J}: 


dR 

dXf 


dU n A ^ dU l d/3 s 
r)\ n + 2^2^ TjflT ' 


( 8 ) 


Solving © yields a globally optimal set of arrival rates for 
the customers. 


B. Local Problems 

We show that the globally optimal solutions yielded by © 
can be found in a distributed way, in which each charging 
statin forms and solves a local problem. In order the lay 
the foundation, we first formulate the local problems in this 
subsection, and relegate establishing the connection between 
the local and global solutions to Section III-C. In order to 
enforce such arrival rates the charging stations compute and 
announce prices for each customer type, which in turn forces 
the customers to operate at the desired (optimal) arrival rates. 
In order to solve this global problem and formalize such 
dynamic between posted prices and adjusted arrival rates, by 
recalling that pj is the prices charged to customers of type j, 
we first define the following local optimization problem for 
each customer n £ {1,..., N } and service type j £ (1,..., J}: 


U n (X n ;(3(X)) 4 U n (X n -f3(X)) - 5>Ar(l-/3 ; ) . 

V 1=1 
gain v_ v _✓ 

cost 

(9) 

By solving this local optimization problem, customer n com¬ 
putes its locally optimal and relevant arrival rates A” = 
[A™,..., A"]. Specifically, the locally optimal arrival rates that 
maximize the local optimization problems {U n } n for all 
Vn £ (1,..., N} and j S {1,..., J} satisfy: 

8U n 

— -Pjii-Pi) =0 . ( 10 ) 

By solving (ITTTb each customer n can compute its own arrival 
rates based on the announced prices. The structure of the 
optimal local arrival rates depends on the selection of the 
utility function U n , which is discussed in more details in 
Section [VI] 


C. Connection Between Local and Global Solutions 


Before proceeding to the details of the relevance of the 
globally optimal solution to these local solutions it is important 
to note that individual customers do not have knowledge about 
the derivatives of their arrival rates with respect to LoLP 
/ 3j which leads to above local optimization problem ll28l . 
Also, the size of the system is fairly large when compared 
to demand of each individual customer, hence each EV does 
not have significant impact on the LoLP. Hence, by comparing 
equations ([8} and (fT0l> we can deduce that by the optimal prices 
for j £ (1,..., J} satisfy: 



N J 

-a-«-‘EE 


l —1 S=1 


dU l df3 s 
d/3 s dXj 


( 11 ) 


When a customer of type j is presented with price p* the 
solution of her local optimization problem © satisfies the 
conditions given in 0, which in turn guarantees that the 
solution of the global optimization problem in ([Sj is equivalent 
to the local ones given in © as established in |28|. Hence, 
the local optimal solutions become the equilibrium since no 
single customer can obtain a higher utility by deviating from 
its locally optimal solution. It is noteworthy that the prices 
{p\,...,p*j} can be considered as congestion prices since 
each customer pays other customers for one unit of marginal 
decrease in their utility due to the rise in the LoLP due to 
increase in the arrival rates. 


IV. Resource Provisioning Framework 

In the resource provisioning problem characterize in © 
the optimal solution C* increases monotonically as the target 
outage upper bounds increase. Hence, one naive approach to 
solve this problem is a linear brute force search, which can be 
computationally prohibitive for large scale system. Instead, we 
aim to provide an establish an analytical connection between 
the optimal solution and the target outage requirements. This 
proposed approach, moreover, is versatile enough to be applied 
to other relevant optimization scenarios, e.g., integration of 
renewable generation or energy storage system and for varying 
customer demand. 

In order to furnish the tools we model the request arrivals 
of the customers as a collection of J independent queues 
(see ll29l ). In this model we start the analysis by first assuming 
that the grid resources are infinite and treat the case of finite 
resources in the next stage. Let Ql x . denote the number of 
customers of type j requesting bj units of power concurrently. 
Further, let S be the sum of offered load by the system that 
is given by 

J 

S = EWc- (12) 

3 =1 

Due to Poisson model for the arrival rates, for the mean and 
the variance of Q ^ we have E [QjJ = var [QjJ = Qj = ^ j- 
Therefore 

j j 

E [S'] = E an( l var [S’] = E (13) 

3 =1 j=1 
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In the case of finite resources C, we can reformulate the 
carried load on the system by defining (}[. as the number of 
simultaneous customers requesting bj units of grid resources. 
Then the blocking (loss-of-load) event occurs if upon an arrival 
of customer j the load on the system is greater than available 
resources Zj=i Qc ’ s greater than C —bj. For a given set of 
charging rates b = [&i,..., bj] and average numbers of arrival 
rates q = [gi,..., qj], let us denote the LoLP by of customer 
type j by /3j : R J —> [0,1], which captures the connection 
between the block probability on one hand hand, and b and 
q , on the other hand. 

By noting that the random number of customers for each 
class are mutually independent Poisson random variables 
we have 


i (q,b) = F\C-b j <'£b j Q 


(14) 


i=i 

J 


= v\c-b j <Y,b j Ql<c\ J2bjQL<c 

i=i i=i 


pj 


bj < ZU b iQ 3 °° 

<c} 

pj 

[Sj=i bjQio < cj 



(15) 

(16) 


In the resource provisioning problem of interest, in order 
to meet multi class QoS targets {Sj}, C* should be at least 
as big as the mean offered load on the system, that is S = 
J2j=i bjQj- However, customers arrive in a stochastic fashion, 
hence it is required to add extra capacity to accommodate the 
fluctuations beyond the average offered load. For this purpose, 
C* is set equal to the mean of the total load E [S'] adjusted by 
an extra term which is a multiple of its variance denoted by x ■ 
var [S], In this formulation more stringent QoS targets lead to 
larger values of x. The objective of the resource provisioning 
problem is to provide a closed form solution to ©. To this 
end, we first scale the system with y > 0 , according to which 
we have the following capacity for the network: 


6 %, x) = ^2 bjqj + X A 
i=i 

and we have the following limiting result ED, fl30): 


A • 

\ 1=1 


lirn v 7 ?/ 3j(q,b ) = 




ip{x) 


(17) 


(18) 


where <p(x) = -^==e x / 2 and ¥>(&) = - 7 = e 4 / 2 dt. 


• 


Now let us define function ib as the inverse of — for Vie, that 

~ C D 


is: 




(19) 

Note that ip(-) is a strictly decreasing function with ip(y)+y > 
0, Vy ( (23), GO) and references therein). Then, the asymptotic 
behavior of the QoS constraint /3 J (q, 6 ) < Sj results: 


x > 


A E b fa 

\ 1=1 


( 20 ) 


The provisioned grid resources should satisfy QoS targets for 
all classes. Hence, the inequality in ( |20) yields 


x > ip 


mm 
1 <j<J 



( 21 ) 


By operating at the lower possible values of x, the minimum 
required amount of resources solving the provisioning problem 
in <[ 6 j is 

EE^l ’ 
1=1 

( 22 ) 

where ip(-) can be easily computed numerically by solv¬ 
ing GD: 


C* = ^b, 

l=i 


Hj 


■ ib | mm , 

1 <j<J bj 



ip(x)J — xVO-bir = 0 , (23) 

and it can be plugged back into (l 22 l) to compute the required 
capacity. 

It is important to notice that providing enough resources to 
meet the QoS targets of the most dominant classes, that are 
the ones with minimum Sj/bj, is sufficient for the remaining 
customer classes. 


-i e -o.5V-(x) 2 _ f _L 

a/2 


V. Computing Loss-of-Load-Probabilities 


Both the pricing-based control problem (summarized in 
section | 6 (c)| with Q}) and the resource provisioning problem 
(summarized in with with ©) fall into multi-dimensional loss 
systems (or multi-rate Erlang loss systems), which are used to 
evaluate the guaranteed performance of networks with limited 
resources. In such systems an arriving customer requesting 
a certain amount of grid resources is either admitted to the 
system or is blocked and we are interested in computing LoLP 
as a function of system parameters. 

Computing the LoLP function requires the analysis of the 
J independent time reversible Markov chains in which the 
state of the system is defined as the number of customers 
of each type, that is Q = [Qc, ,Qc\ an d the state space 
is denoted by H = {Q : Y^j =1 bjQio < C}. Let Q :i c denote 
the maximum number of customers of type j that can be 
served simultaneously. Assuming the order &i > ... > bj > 0, 
without loss of generality, provides that 0 < Qq < ... < Qc- 
Then the probability of being at state Q is (29]: 


r(Q) = n 


o-1j 


7 = 1 O' 


•So! 


(24) 


Next, similar to (fl5l > we condition on the finite capacity, 
and compute a generic state Q probability distribution as: 


7t(Q) 


?t(Q) 

J2qc n ^(Q) 


(25) 


Next, let us define the blocking states (LoLP) for customer 
type j as 


^ J ={Q:C-b j < J2 b kQc < C} ■ 

k =1 
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Hence, ( fl6t can be re-written as: 

^( q , b ) = Y tt(s) = !- Y 7r ( s ) > (26) 

setj sg'S’j 


where the second term represents the probability that the 
charging station has a total capacity of C—bj instead of C and 
7 r(s) is the steady state probability mass function. Furthermore, 
let us define function H(C, J) as 


H{C, J) 


= e n 

{Q: bQ < C} 



3= i Q j '- ’ 


(27) 


based on which one can compute the LoLP for class j 
explicitly as: 


Pj(q-b) = i - 


H(C-bj,J) 
H(C, J) 


(28) 


The set {Q : bQ < C} in (f27l) contains all the states corre¬ 
sponding to which no outage occurs. While (l28l > provides an 
explicit representation for j3j, the associated computation can 
be costly when the system capacity C is large. Considering a 
real-world scenario in which the number of classes J is typi¬ 
cally varies between 3 and 5, and C is in the order of Mega¬ 
watts, computing LoLP could be easily carried out via the 
Kaufman-Roberts algorithm (Algorithm [TJ which involves a 
simple recursion which computes the occupied resources ll24l . 
Let c denote the amount of resources being in use and 


a(c) = P {c units of power in use} . 


We subsequently have 


*(c) = 


1 


{Q: bQ < c} 


Q ] H(C, J) 


(29) 


and the LoLP corresponding to customer type j can be 
calculated by using 

bj~ i 

Pj (q, b )=Y a ( C ~ *) • ( 30 ) 

2—0 

Note that the above derivations are based on the assumption 
that grid resources are discretized (e.g., 1 kW is considered 
as 1000 discrete serving units) and the interpretation of the 
Algorithm [T] is that LoLP of customer type j equals to 
some of occupied states that are in the range of C — (b :l — 
1C — 1,C. Another important measure of interest for 
solving the optimal pricing problem given in ( fill is the set 
of derivatives of LoLP with respect to offered load by each 
customer class. In the proposed multi-class customer model, 
there is no explicit formulae for performance measures (LoLP) 
in terms of input parameters (C, A j,Pj, etc.). In this paper we 
follow the methods based on convolution algorithms Il32l to 
compute the derivatives of the LoLP. To that end the derivative 
of the LoLP associated with customer type j with respect to 
the traffic intensity of another class j\ ^ j 2 can be computed 
by using the function a(-) as follows. 

^ = a(C - b h ) + a(C - b h - 1) + • • • (31) 

dq j2 

+a(C - b h - b h - 1) - (1 - ft*) ft' 


Algorithm 1 Kaufman-Roberts Algorithm lf24l 

Set k( 0) = 0 and n(i) = 0 for i 6 1R~ 

for i=l to C do 

= tE ~bj) 

end for 

Compute H = X^iLi K (*) 

for i =0 to J do 

«(0 = # 

end for 

for j =1 to J do 

Pj(q>b) = Ei^c-bj+i a (*) 

end for 


A useful property of the LoLP is the elasticity property, which 
entails that the sensitivity of customer LoLP j to customer 
class k is the same as the sensitivity of customer LoLP k to 
type j lf33ll . This property is given as yY = ^ 2 .. 

Notice that Kaufman-Roberts algorithm is used in both of 
the proposed frameworks and the complexity of this algorithm 
is O (CJ) which provides a considerable improvement over 
solving the LoLP through (l28l) which has a complexity of 
O (C J ) ll24l . l25l . For the pricing-based control problem, 
we compute the derivatives of multi-rate blocking probabil¬ 
ities using the convolution-based algorithm [32] and similar 
to Kaufman-Roberts algorithm, the order of complexity is 
O (CJ). 

VI. Simulation Results 
A. Case Study-I: Pricing-Based Control 

In this subsection we present a case study to evaluate 
the social welfare maximization problem of the pricing-based 
control framework presented in section [III] The parameter 
setting for our case study is as follows. For the ease of 
representation we assume that there are two customer types, 
namely, fast charging customers (type-I) and slow charging 
customers (type-II). In order to mimic the typical charging 
rate for a fast DC charging (that is 50kW), we set &i=50 
units. Also, charging duration takes around 20min, so mean 
service rate is set as /ii=3. In a similar manner we tune the 
parameters for the slow charging customers as b 2 =7 units and 
/i2=0.42. For the utility functions we adopt the widely-used 
logarithmic utility function ll 22 l . fST\ . The utility function of 
a single customer increases with the arrival rate and decreases 
with the customer LoLP according to 

f ZU “j log(l + Xj) - 9j log(l + P 3 C ) A j > 0 

H 

to Xj < 0 

(32) 

where ujj and 6j are the weights of each class. Note that 
the weight is higher for customer types with higher demand 
bj, mainly because fast charging customers demand more 
resources and hence gains more utility. A simple example is 
presented to clarify the matters. Assume that system capacity 
C = 500, and the weights are chose as wi=20, W 2=10 and 
0i=6O, 02 = 20 . Then the social welfare could be maximized by 
setting the arrival rates to Ai=8.6638 and A 2 = 5 . 2001 . For this 
arrival rates prices pi and p 2 are computed to be 0.3197 and 
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(a) Social welfare computation. 




(b) Optimal arrival rates (Ai^) with two (c) Optimal arrival rates (Ai, 2 ) with two 
classes. classes. 


Fig. 6: Numerical Evaluation-I: system capacity C i6 fc=450 + 50 sin(27rfc/80), k G (0,10,..., 80) and C is assumed constant 
for a duration of T=10 


0.2211 and the resulting LoLPs are /3 1 =0.0097 and /1 3 =0.0009, 
and the maximum utility is computed to be 59.1238 units. 
Notice that u>j is chosen greater than 0 :t so that station operator 
is motivated to provide a good level of QoS. We further explore 
the relationship between the optimal arrival rates and the social 
welfare. As shown in Fig. |6(a)| if the arrival rates deviate from 
its optimal value, the social welfare reduces. 

We proceed to provide more numerical evaluations. In 
the first setting we present the relationship between sys¬ 
tem capacity and optimal arrival rates for time-dependent 
case, where station capacity varies over time and follows 
C(fc)=450 + 50 sin(27r/c/80) (due to grid conditions) and 
also C{k ) is assumed to be constant for every T = 10 
duration. Results depicted in Fig. |6(b)| shows that, it is more 
beneficial to accept fast charging customers as they improve 
the social welfare function more than slow charging ones. 
For the given set of arrival rates, the corresponding prices 
in (fTTb and loss-of-load-probabilities are given in Figs. |6(c)| 
and |7(a)| respectively. Obviously, since fast charging uses more 
resources, the corresponding prices are higher and due to 
higher arrival rate the corresponding LoLP is higher than the 
slow charging customers. Finally, we present the correspond¬ 
ing social welfare in Fig. |7(b)| 

B. Case Study-II: Capacity Planning 

We proceed to compute minimum amount of grid resources 
to provide QoS guarantees to two customer classes with the 
same set of parameters (pi=3, ^2=0-42, 6i=50, &2=7) for a 
wide range of QoS targets (0.001 < <5i < 0.05, 0.001 < 62 < 
0.05) with fixed arrival rate Ai i 2 =5. The results depicted in 
Fig. |7(c)| presents that, since type-1 is the dominant class in 
(most) regions where 5\/b\ < 82 / 82 , providing resources for 
class-1 already satisfies the QoS targets for class-2. Next we 
compute the required capacity for different range of arrival 
rates and fixed QoS targets <5i=^2=0.03. Results depicted in 
Fig. |8(a)| can be used as a guideline to choose the required 
capacity for given arrival rates. 

We proceed to investigate the percentage of reduction in 
station capacity (C) for a given LoLP targets. The motivation 
is that instead of providing zero percent LoLP, station operator 


can sacrifice to reject small amount of customers and reduce 
the stress on the grid. We use the same parameter setting as 
above and computed the required capacity respect to almost 
zero LoLP (<5 i.2 =10~ 6 ). As presented in Pig. |8(b)| even one 
percent QoS targets leads to significant savings in station 
capacity. Our final evaluation is on evaluating the system 
performance for non-homogenous arrival rates. Let’s assume 
that arrival rate for class-II is constant and A 2 (fc )=10 and 
Ai=10 + 2 sin(27rfc/80). The QoS targets are set as (5i=0.04 
and 52=0.01. Then station operator can provision the system 
according to peak hour demand (computed as (7=683) and 
hence meet the QoS requirements at all times. Results are 
presented in Pig. |8(c)| 

VII. Concluding Remarks 

In this paper we provided two important design problems 
for electric vehicle charging stations with multiple classes of 
customers. We considered two cases. In the first one, given a 
capacity of a station and infinite amount of customer request, 
we provide a framework to compute the optimal arrival rates 
such that the total social welfare is maximized. This is a very 
typical case for charging stations located in big cities. On the 
other hand, in the second case we considered charging stations 
located in small cities. This time our primary concern was to 
calculate the minimum amount of grid resources such that each 
customer class is ensured with a certain level of QoS targets 
(LoLP). 

This initial work can be expanded in different directions. Lor 
example each station can employ an energy storage system that 
can aid to further reduce the strain on the power grid. The ESS 
can be charged during light traffic and the stored energy can 
be used to meet EV demand during peak hours. Another future 
research interest would be to consider a different resource 
policy to optimize the resource usage (e.g., portioning for 
each class etc.). In this paper we assume one micro-grid level 
charging. Another research direction can include a customer 
routing among the micro grids or stations. 

In this work, we have concentrated our focus on a network 
of stations fed by a single substation. To that end, our final 
research direction is on considering a general network case 
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(a) LOLP performance. 



(b) Social welfare computation. 

Fig. 7: Numerical Evaluation-II 



QoS Target 0.05 0 .05 004 QqS Target g 


(c) Capacity planning for varying QoS targets 
with two classes. 



(a) Resource provisioning (fixed Si = 82 = 
0.03) 



Fig. 8: Numerical Evaluation III 



Time 

(c) LoLP Performance, Cmin = 683. 


and address the congestion issues in a grid composed of power 
lines and other elements with different capacity ratings. 
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